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PREFACE 


Part  of  the  research  program  of  The  RAND  Corporation 


consists  of  basic  supporting  studies  In  mathematics.  One 


aspect  of  this  is  concerned  with  the  solution  of  partial 


differential  equations.  In  this  field,  the  technique  of 


difference  approximations  has  been  extremely  useful. 


In  the  present  Memorandum,  the  applicability  of 


higher-order  difference  approximations  of  an  unconventional 


type  to  the  solution  of  partial  differential  equations  Is 


Investigated. 
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In  this  Memorandum  techniques  of  approximating 
the  solution  to  a  partial  differential  equation  are  in¬ 
vestigated.  Using  the  first  of  these,  the  solution  is 
approximated  at  each  stage  by  a  higher-order  difference 
algorithm.  The  second  technique  is  that  of  storing  the 
function  at  each  stage  by  a  polynomial. 

Numerical  results  were  obtained  from  a  FORTRAN  pro¬ 
gram  applying  these  techniques  to  the  equation  u^  ■  uu  . 
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HIGHER-ORDER  APPROXIMATIONS  TO  THE  COMPUTATIONAL 
SOLUTION  OF  PARTIAL  DIFFERENTIAL  EQUATIONS 

1.  INTRODUCTION 

It  has  been  shown  In  [1]  that  certain  partial  dl£> 
ferentlal  equations  could  be  solved  by  an  approach  which 
has  some  advantages  over  the  usual  methods  of  approxlma" 
tlon  by  difference  equations.  This  approach  was  shown 
to  be  superior  In  the  following  ways:  1)  the  solution 
Is  re-created  at  desired  points  by  an  approximating  poly¬ 
nomial;  and  2)  the  algorithm  for  approximating  the  solu¬ 
tion  exhibits  certain  desired  properties  of  the  actual 
solution.  To  Illustrate  this  method  It  was  shown  In  [2] 
that  the  equation 

u^  -  uu^  ,  u(x,0)  -  g(x)  ,  (1.1) 

which  possesses  the  analytic  solution 

u  -  g(x  +  ut)  ,  (1.2) 

where  g(x)  Is  here  assumed  to  be  an  odd  function  of  period 
two,  can  be  approximated  by  the  algorithm 
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u(x,t  +  A)  -  u(x  +  u(x,t)A,t)  . 


The  variable  t  was  constrained  to  values  t  «  0,A,2A,..., 
but  X  assumed  arbitrary  values  in  the  interval  [0,13.  The 
algorithm  (1.3)  clearly  preserves  boundedness  and  non¬ 


negativity,  and  upon  expanding  in  a  Taylor  series  it  is 

.2 


seen  that  the  approximation  is  accurate  to  0(A  ).  In 
[2]  it  was  suggested  that  more  accurate  results  could  be 
obtained  using  a  higher-order  algorithm.  It  is  the  pur¬ 
pose  of  this  paper  to  investigate  the  effect  of  improving 


the  approximation  to  (1.1)  by  using  an  algorithm  accurate 


to  0(A  ). 


2.  HIGHER-ORDER  APPROXIMATION  , 
Let  the  algorithm  be  given  by 


u(x,t  +  A)  ■  u(x  +  u(x  +  u(x,t)A ,t)A ,t) 


(2.1) 


which  also  preserves  boundedness  and  non-negativity.  Ex- 

,3, 


panding  both  sides  of  (2.1)  in  a  Taylor  series  to  0(A  ) 
and  applying  the  relations 


u^^  ■  uu  +  uu  ^ 
tt  X  xt 


(2.2) 


A  ^  ^ 

+  ~  +  O(A^)  -  (2uu^^  +  u^u^)  +  O(A^) 


holds,  and  hence,  (2.1)  is  accurate  to  O(A^), 

Again,  let  t  ■  0,^,2 A,...,  and  at  each  stage  of  the 
calculation  let  u(x,t)  be  obtained  by  means  of  the  finite 
sum 

M 

u(x,t)  a:  ^  Ujj(t)  sin  (nirx)  ,  (2.4) 

n-1 


(2.6) 
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H6nc6|  va1u68  u{k/Rft) •  k  ■  1|2, •  •  •  store  u(X|t) 

at  time  t,  and  by  way  of  (2.1)  u(x,t+A)  can  be  obtained. 


3.  NUMERICAL  RESULTS 

A  FORTRAN  program  for  the  IBM  7090  was  written  to 
compare  results  obtained  using  (2.1)  with  those  using  (1.3). 

(a)  g(x)  -  0.1  sin  ffx,  0  ^  x  s  1 
where  M  -  R  -  10,  A  -  0.1,  0  ^  t  ^  10 

yyition  requires  22  seconds  using  (1.3)  and  27 


seconds  using  (2.TT?**‘*^Afi%5j[;^own  in  [2],  using  (1.3)  gives 
errors  in  atlc,.t)  ^  .003  for  times  t  Tn“''l"’*4iiiyii|||,Miii  1 1  As 
shown  in  the  following  table,  using  (2.1)  usually  gives'^* 
slightly  better  results  than  (1.3)  for  times  t  less  than 
three . 


X 

t 

exact 

<1.3) 

error 

(2.1) 

error 

.1 

1 

.043598 

.043058 

.00054 

.043580 

.00018 

.3 

2 

.lOOOOO 

.100097 

.00010 

.100159 

.00016 

.7 

2 

.055795 

.055416 

.00038 

.055854 

.00006 

.9 

1 

.023732 

.023579 

.00015 

.023736 

.00000 

For  t  ~  10 /ir,  the  time  at  which  the  shock  occurs,  the 
values  of  u(x,t)  using  (1.3)  were  in  error  about  3  per 
cent  at  X  2  0.2,  0.4,  0.6,  and  0.8.  Also,  the  largest  °  'Tl 
errors  are  made  for  small  x,  since  the  shock  occurs  at  H 


X  “  0.  It  is  surprising  to  note  that  (2.1)  gives  poorer 
results  than  (1.3)  for  time  t  -  3.1  (see  following  table). 


X 

t 

'  exact 

(1.3) 

error 

(2.1) 

error 

.1 

3.1 

.094349 

.099887 

.00554 

.101926 

.00758 

.7 

3.1 

.046784 

.047194 

.00041 

.047843 

.00106 

.9 

3.0 

.016137 

.016212 

.00008 

.016456 

.00032 

(b)  g(x)  -  0.1  sin  ffx,  0  s  X  s  1 
where  M  ■  R  »  10,  :  A  -  0.05,  0  ^  t  ^  10 

To  determine  the  reason  that  a  higher-order  algorithm 
.do|W,^not  give  bet^^  near  the  ^ock,  the  same  cal- 

culation  was  done  using  a  smaller  stepsize  A.  For  times 
t  less  than  three,  the  accuracy  in  u  calculated  using  (1.3) 
is  generally  better  for  a  smaller  A;  however,  the  calcu¬ 
lation  using  (2.1)  is  inconsistent:  sometimes  better,  other 
times  worse.  Finally,  using  (2.1)  and  a  smaller  stepsize 
does  not  give  better  results  for  t  ~  10/ff. 


.  'S..,,,,  tn,,;question  the  value  x. . 

u(x+u(x  ,.t)A, t)  is  the  polynomial 

given  by  (2.4).  However,  this  polynomial  .i's  used  a  sec¬ 
ond  time  in  calculating  -x.^  -  u(x,t+A)  «  u(x+Xj^A,t)  in  the 
case  of  the  higher-order  algorithm.  Now,  if  tht^  poly-’ . .. 
nomial  approximation  error  is  sufficiently  large,  then 


lncrMt«d.  This  is  ths  reason  for  the  poor  results  near 


the  shock  in  examples  (a)  and  (b),  above. 

(c)  g(x)  ■  0.1  sin  wxf  0  ^  X  ^  1 
where  M  ■  R  ■  20,  ^  ■  0.05,  0  ^  t  ^  10 


exact 


error 


3.1  .094349 

.7  I  3.1  .046784 


Execution 

Time 


a. 3) 


.043320 

.099997 

.055591 

.024228 

.091836 


.046269  n  .00051 


87  sec 


r2.n 


error 


.043593  .00000 
.099995  .00000 
.055795  .00000 


.023733 

.091602 


.00000 

.00275 


.046397  I  .00038 


123  sec 


In  example  (c),  however,  a  20^  degree  polynomial  is 

used.  As  seen  in  the  above  table,  better  results  are  now 

3 

obtained.  For  times  less  than  three,  the  0(A  )  algorithm 
gives  results  precisely  to  the  flftli  place,  while  for 
t  ■  lO/s  large  errors  (mainly  those  near  se  «  0)  are  reduced. 


4.  CONCLUSION 


— ‘  It  appears  from  the  results  that  in  order  to  reduce 

:  ^e._.computational  error  by  increaslt^t  t^e  order  of  the 
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